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Abstract — In this paper, an algorithm is presented for 
identification of real eigenvalues on right half of the s-plane, 
for linear systems, hence determining instability of the system. 
The proposed approach is based on Gerschgorin theorem and 
a new approach of Bisection method. The method is efficient 
since there is no need to determine all real eigenvalues and 
also characteristic polynomial of the system matrix. It has been 
found that in some class of control system problems, the method 
needs minor computations. The proposed approach is useful, 
particularly, in power system applications, where the order of 
system is large. A power system problem and numerical 
examples are illustrated using proposed algorithm. 
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I. Introduction 

The stability problems arise mainly in the power 
system whenever perturbation occurs. These perturbations may 
occur because of change in the parameters of the system. The 
system response following the perturbation may be either stable 
or unstable. The control engineers are always interested in 
whether or not the system is stable. In practice, the systems are 
nonlinear in nature. The analysis and synthesis of nonlinear 
system are quiet difficult. So, these nonlinear systems are to be 
linearised around the operating point to obtain linear state 
variable model described by the following state variable 
equation, 

#= Ax+ Bu ■ (1) 

Where, 

A = System matrix, g = Input matrix, 
X = State vector, u = Control vector. 

System matrix A consists of parameters of the 
system and its eigenvalues play very important role in the 
stability of the system, particularly in the power system. 
From eq. (1), the eigenvalues of the matrix A gives the system 
behavior, whether the system is stable or unstable and if 
stable, how much it is relatively stable. Thus, the eigenvalues 
are in general, functions of all control and design parameters. 
The stability of the system can be determined by applying 
Routh's criterion to the characteristic polynomial of the 
system [1]. This criterion gives the presence of eigenvalues 
of the system on the right-half of the s-plane. In the similar 
manner, the proposed method tries to determine the presence 



of real eigenvalues of the system matrix A belonging to the 
right half of the s-plane without computing the actual 
characteristics polynomial and eigenvalues. In [2] , technique is 
presented to identify real eigenvalues using Gerschgorin 
theorem [3]. But, the algorithm presented in [2] fails when the 
eigenvalues of the system are of repetitive nature and therefore 
the system which is actually unstable, may show stable by the 
existing algorithm [2]. Today, fast computing software such as 
Matlab is available by which stability of any system can be 
tested very easily. But, in some cases, because of rounding of 
errors, truncation errors and ill conditioning, results shown by 
Matlab may be inaccurate. Hence, in such class of problems, 
Gerschgorin bounds will be helpful, because based on the 
location of the bounds, we can decide the stability of system. 
Moreover, the calculations of bounds require very minor 
computations in comparison with the calculation of actual 
eigenvalues. 

Thus, in this paper, an algorithm is presented to check 
the instability of the system matrix using Gerschgorin theorem 
[5-7] and a new approach of Bisection method [8]. In order to 
show the effectiveness of the proposed method, the same power 
system model [4] is considered which has been considered in 
[2] and computational efficiency is highlighted. Similarly, 
various other examples are also illustrated. 

In this paper, the following notations are used in 
mathematical developments. Complex plane is dented by £ ; 

Open left-half plane is denoted by £ - ; open right-half plane 

is denoted by £ + ; belong to is denoted as I . 

II. GERSCHGORIN THEOREM 

For 77 J 2, let A be any 77 77 real matrix (written 
A= [a,. y ]i j rin ,i,j = 1,2 »), and let s (A) 

denotes its spectrum Q^^^X^LL^LLdLl^LijJ^jQ^; 
A familiar result of Gerschgorin as follows. 

Theorem II. 1: 

The largest eigenvalue in modulus of square matrix 
A cannot exceed the largest sum of the module of the elements 
along any row or any column. 
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Since the eigenvalues of j^ T are same as those of A , therefore, 

the theorem is also true for columns. 
Theorem II. 2 : 

Let p k be the sum of moduli of the elements along 

the fc tH row excluding the diagonal element a k k . Then every 

eigenvalue of lies inside or on the boundary of at least one of 
the circles in s-plane. 
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p t ,k= l,2,...,n, (3) 
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where, fl, , is center of circle with radius p t . By 
denoting each of these 77 circles by R k , we get, 

/? t =|/-%|}=^=U n (4) 

Therefore, each of the eigenvalues of matrix A must lie 
in the union of ' S ' of these W circles in the s-plane, i.e^ 



s=m 



(5) 



*-l 



It is also noted that A and A have the same 

eigenvalues, Thus by applying theorem U.2 to A , it 

u 
yields another set S . Thus, 

S(A)= SI S (6) 

In case of real symmetric matrix A and A are same. 
Thus, 

s(A)=So L S. (7) 

The eigenvalues of matrix A will lie in the intersection 

of these two unions S and S in the s-plane. Thus, by 
considering the intersection of these Gerschgorin circles, 

bounds under which real eigenvalues will fall, are obtained. 
The bounds obtained are nothing but the extreme ends of the 
intersection of Gerschgorin circles. Let these bounds are denoted 
by E and D, where D is the extreme left bound and E is the 
extreme right bound [6-7]. 

If D, EI £ , then "T I £ ,i= 1,2,..., n, (all the 
eigenvalues lying on the left of the s-plane), system is stable. 



If D, EI £ + , then"/, I £ + ,/= 1,2,. ..,n, (all the 

eigenvalues lying on the right half of the s-plane), system is 
unstable. 

Remark 1: The above theorem is also true for the 
complex elements in the system matrix. 

m. MAIN RESULTS 
Theorem III. I: 

When there exist odd number of eigenvalues on £ + 

for a given system matrix A j ■ " " , then, 



\ll- A[__ <0. 



(8) 



Proof: Consider A as any n' n real matrix 
(written A=[fl. ]I j ™,i,j= l,2,...,n), and 
/■I £,/= 1,2,... ,71 be the eigenvalues of the matrix. Now, 

\ll- A\= /(/)=(/- /,)(/- l 2 )...(l- l„).(9) 
From eq. (9), depending upon the value of/ , we get, 
/(/ ) / = () > ,for"/ ; I £,£ c( ' +) /= 1,2,. ..,71, (10) 
where, £ e (" •+ ) indicates only complex conjugate eigenvalues 
on £ ■ and £ + . 

f0) l=0 > 0,for"/,i £-,, ,1= 1,2,. ..,71 ,(11) 

and for " / ; I £ + , j , i = 1, 2, ..., n , we get, 

/(/ ) l=0 > jf /= 2,4, 6,..., n, i.e., n 

and 



is even 
(12) 



fQ ) 



1=0 



if /= 1,3,5,..., 77, i.e., 77 is odd. 

(13) 
Thus, from above eq. (9-13), it is proved that when there 
exist odd number of real eigenvalues on £ , then. 



I"- 4=o <0 



(14) 



Theorem III. 2: 

Suppose for a given matrix,, 

At i '\D I =/(/} o0 =M-'4- md 

D,= /(/ ) + = / /- A\, + , then there exist at 

L ' t =e I \l=e 

least one eigenvalues G £ , when Aj > A 2 . 

Proof: Consider a continuous function / (/ ) which include 
all the real and complex eigenvalues, lying on £ and^ 
£ . It can be written as 
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Where, 



/(/„.)=(/- xjl - x j+l )...(l - x p ), (16) 
/(/^)=(/-a, + 1 )(/- V2 )...(/-.:)J17) 



Ov +3 ) 2 )...((/-^) 2 + (v s ) 2 X 

f(l r ,)=(il-x s J + (y s J)((!-x s J + 



()^fU(l-x n f + (y n f). 



(18) 



(19) 



/(/ ), /(/ _ ), shows the product of all the real 



nr 



eigenvalues on Jh£^ Similarly, 

/'(/ ^ ) and /'(/ ) shows the product of all complex 



cC 



conjugate eigenvalues on the £ " and £ respectively. 
Therefore, from eq.(15) and eq.(16-19), the following 
conclusion can be drawn. 

if/I i.£-,then/(/ rC .^ e+ >/(/ rC .) i=0 , (20) 



if /I £,£-, then f(/ r \ + > f(l _),_ , (21) 



jf.il i ,£ + ,then/(/ .)_,</(/ A , (22) 



rC 1 



rC + 'l- 



if/i£,£ + ,then/(/ cC+ ) ; = e+ </(/ cC+ ) i = .(23) 

From above eq. (20-23), when /(/ ) ;=0 > /(/ ) i=g , 

then there will be existence of eigenvalues either real or 

complex, on the £ , but the converse may not be true. 
This completes the proof. 

Theorem HI. 3 : 

For a continuous function f(l), if there is a 

repeated (double) root at t d , d and d be a small value in 

its neighborhood on both sides and assuming that there is no 
eigenvalue between d x and d y except double root at t d 

\{{t d - d 1 )>{t d - ^)and/(/) !=vdj >0,then, 
/(O J -^ > 0md/(/),. (jf . dj >/(/) N ^(24) 

Proof: Suppose, there is a real ejgerrralue t I £ , d r 
Mid d be a small value in its neighborhood, then using 
conventional Bisection technique, 



jf/(/ r 4)>0 then f(t l +d y )<0, 



if /(/,- d x )<0 then f(t l +d y )>0. 



(25) 



(26) 
But, when there is a repeated root at t d and by 

assuming small values d x and d in its neighborhood, the 

above equations are not satisfied. It is because there are two 
sign changes at the same moment, therefore, effectively 
there is no sign change (ex. For first root, - sign change to + 
and for the second root, - sign change to +, so, finally, two 
negative sign changes to positive sign and two positive sign 
change to positive sign). But, the value of the determinant 
decreases while approaching towards double root and its 
value becomes zero at double root and there after the 
value of determinant increase or decrease either in positive 
or negative direction, having same sign which is before 

crossing the root. Thus, for double root at_T d , if 
(V 4)> i f d- ^ X then, 



and 



f{r d - 4)>0,/(* d + 4)>0 (27) 

/('Wd >f(!),- tj+ <L- (28) 



This proves the theorem. 

Remark2: We have only considered the case of repeated 
eigenvalue when/^- (^.)>0.1t is because when 

f(j d - cL ) < 0, then the system is unstable. 

Remark 3: While moving / from to E in suitable 
number of steps, whenever / / - A\ is negative or there 

is a decrease in its value at any instant, then the system is 
always unstable. Based on the concepts of Gerschgorin 
bounds and using above theorems, we are presenting 

algorithm for identification of real eigenvalues on the £ , 
hence the instability of the system. 

IV. ALGORITHM 

Using above theorems and Gerschgorin bounds, a step 
by step procedure for determining instability of the system 
matrix ^4 is as given below. 

Step 1 : Enter order n and elements of the 

matrix, A = [a ■ • ] I j 

Step 2: Calculate bounds D and E using Gerschgorin 
Theorem. 

Step 3: If D, E < 0, then the system is stable, 
otherwise, go to next step. 

Step 4: If D, E > 0, then the system is unstable, otherwise, 
go to next step. 

Step 5: IfD<0andE = 0, thenifD = |/ /- A[ =Q = 0, 
then system is unstable, otherwise, the system is stable. 
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If D < and E > 0, or D > and E < 0, then the following steps 
are applicable. 

Step6: For/ = Q, we calculate, D = / /- A ,_„, 

if D < , then the system is unstable, otherwise, go to next 
step. 

Step7:For/ = e + , we calculate, D t = Z /- A \,_ + , 

if D j < or D > p j , then it indicates the existence 

of at least one eigenvalue belong to £ + . Hence, the system 
is unstable, otherwise, go to next step. 

Step 8: If there are sign changes in the values of 

D= //- A when X is varying from zero to E in H number 
of steps, or there is a decrease in the value of the determinant 
D at any instant, or there is repeated eigenvalue as per the 
theorem III. 3 ,then, it indicates the existence of real eigenvalues 
belong to . Thus, the system is unstable. 

V. APPLICATIONS 
V.l Example: 
Consider the matrix A as [7] 



[A]~- 



Now applying Gerschgorin theorem to above matrix , 
we calculate bounds, as D=-12 and E=-l. Since the bounds are 
completely lying on the left half of the s-plane, the system is 
stable. The Gerschgorin circles and bounds are shown in fig. 1 . 
The actual eigenvalues of the system matrix are -3.0000, - 
7.0000-3.3944, -10.6056. The stability has been decided without 
computing the characteristic polynomial and eigenvalues, hence 
require minor computations. For such class of problems, 
proposed analysis is efficient in comparison with the Matlab 
where in eig (a) command is used to identify eigenvalues and 
hence the stability. 
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Fig.l. Gerschgorin circles and bounds for example V.l 

V. 2 Example 

Consider the system matrix [2, 4] 

"-1.18 0.985 






314 

-0.341 -0.5 0.375 

48.1 -107.5 -205 -54 

295 295 -170 -638 



— I 



0.1125 

431.5 

-19.65 



Applying Gerschgorin theorem to above matrix A > we 
get D=-680 and E=680, as shown in fig. 2. Now, we apply 
proposed algorithm to the above matrix. Using step 5 of the 
algorithm, we get,D D = I' /_ A \ i-o = - 6116382. Since, the 
determinant is negative, the system is unstable. The actual 
eigenvalues of the matrix A are -.7.52,- 0.35,- 8.5,- 37.04± /25.The 
proposed analysis needs only one iteration, whereas as the 
existing method by pusadkar et al. [2] needs larger number 
of iterations. 



t5C0 r 




■1000 -500 500 1000 1500 

RmIMk> 

Fig. 2. Gerschgorin circles and bounds for example V.2. 
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V. 3 Example: (Identification of repeated eigenvalues) 
Consider the matrix A as 

"0 1 

1 



[-■!] = 



-25 



■15 



The Gerschgorin circles and bounds for the above 
matrix A are E =25 and D=-25 as shown in fig. 3. 




Fig. 3 . Gerschgorin circles for Matrix [A] 

Now we check the stability using the proposed 
algorithm. The / /- A for various values of / are as given 
below in table 1 . 

Table 1: Identification of repeated eigenvalues 



Steps 


/ 


\1 1- A\ 


1 





25.00 


2 


0.1 


26.4110 


3 


2.08 


26.22 


4 


4.16 


3.58 


5 


6.24 


11.13 


6 


8.32 


102.72 



It is observed from Table 1, that the value of determinant 
decreases from step 2 to step3. Hence, the system is unstable, 
and there is no need for further iterations to check stability. 



But, in order to check the repeated eigenvalue, we have 
shown further iterations. From step 4 to step 5 in the Table 1, 
instead of decrease in the value of determinant, there is an 
increase in its value. So, from theorem III. 3, the value of 
determinant is increased after a decrease, without change in 
sign, hence it is concluded that there exist repeated eigenvalue 
between 4. 16 and 6.24 on the real axis of the s-plane. This shows 
that the system is unstable. The approximate repeated eigenvalue 
is the average of 4. 16 and 6.24, i.e., 5.2. The actual eigenvalues 
of the above matrix are -1,5, and 5. Thus, the drawback of the 
paper [2] is improved in the proposed analysis. 

VI. CONCLUSION 

The proposed method is useful for determining the 
instability of the system when there exist real eigenvalues on 
the right half of the s-plane. It is based on the Gerschgorin 
theorem and a new approach of Bisection method. The method 
is extremely efficient when bounds E and D lie completely on 
either side of the s-plane or when there exist odd number of 

real eigenvalues lying on £ + . In future, an algorithm can be 
developed for the identification of the complex conjugate 
eigenvalues which lie on . This new method will be used as an 
alternative to Routh Hurwitz criteria for the stability of the 
system. 
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